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INTitODUCTION 

The  classical  ohermodynanic  iunctions  of  a  system  are  determined 
from  th3  observable,  macroscopic  properties  of  a  system.  The  macro- 
scopic properties  in  turn  may  also  be  determined  from  the  micro- 
scopic point  of  view  by  considering  the  properties  and  statistical 
behavior  of  the  atoms  and/or  molecules  in  the  system.  The  knowledge 
of  molecules  and  the  statistical  characteristics  necessary  for  the 
evaluation  of  the  thermodynamic  functions  is  incorporated  in  the 
partition  function,  and  an  underatandia^^  of  the  partition  function 
permits  the  derivation  of  all  the  thermodynamic  properties  of  a 
gaseous  system  in  the  statistical  approach. 

The  partition  function,  by  definition,  is 

£i 


»« 


KT 


Q  =  Is^e        -■'■  (1) 

1=0  1 

vrhere  Q  is  the  partition  function,  g.  is  the  statistical  weight,  ^  . 
is  a  quaiixum-mechanically  allowed  energy  level  of  the  molecules  which 
constitute  the  system,  K  is  the  Boltzmann  constant,  and  T  is  the 
absolute  temperature  of  the  system. 

Physically,  the  partition  function  is  a  pure  number  and  might 
be  evaluated  for  a  system  in  any  physical  state,  containing  any 
number  of  molecules.  It  is  a  well  defined  function  used  to  derive 
thermodynamic  properties  of  a  system  from  the  state  of  atoms  or  mole- 
cules. As  an  example,  the  five  basic  thermodynamic  properties  S,  E, 
PI,  A,  and  G,  and  the  two  derived  thermodynamic  properties,  Cy  and  C  , 
of  a  gaseous  system  can  all  be  expressed  in  term.s  of  the  partition 
function.  The  definition  ox'  the  partition  function  is  the  result  of 
natural  grouping  of  variables  in  the  derivation  of  the  Maxwell-Boltz- 
mann  distribution  law,  the  law  which  describes  the  statistical  be- 
havior of  any  molecular  gaseous  system.  It  can  be  seen  from  the 
following  sequence  in  this  report  that  all  the  thermodynamic  proper- 
ties of  a  molecular  gaseous  system  can  be  calculated  on  the  basis 


of  a  Kaxwellian  model  ".-.'•ith  knovm  speotroscopio  information.  The  value 
of  this  kind  of  calculation  lies  in  tvro  facts;  first,  the  thermodynamic 
properties  of  a  system  that  may  be  unobtainable  by  macroscopic  experi- 
ment can  be  evaluated  using  this  method,  and  second,  even  though  the 
macroscopic  experiment  can  be  carried  out,  the  statistical  method 
yields  thermodyna'nxc  information  with  a  higher  accuracy  and  precision 
thcin  the  classical  method. 


I.  Energy  Level  6 ,  and  Statistical 
Weight  g,  in  a  System 


Theoretical  aiid  experimental  investigations,  from  the  areas 
of  quantum  mechanics  and  spectroscopy,  have  yielded  the  pictures  of 
the  structure  of  individual  molecules  in  the  ideal  gaseous  state, 
that  is,  for  the  molecules  at  pressures  sufficiently  low  that  inter- 
actions between  the  molecules  are  negligible.  From  the  Schrodinger 
wave  equation  (Ref.  6), 

H  H'  =  E  4^  ^^^ 

where  H  is  the  Hamiltonian  operator,  ^  is  the  wave  function,  i.e.  the 
eigenfunction,  and  E..  is  the  energy,  i.e.  the  eigenvalue,  the  well 
defined  energy  levels  are  established  theoretically.  Experimently, 
the  energy  of  a  system  can  be  found  from  the  observed  spectroscopic 
data.  Assuming  that  there  is  no  interaction  between  the  translational, 
rotational,  vibrational,  and  electronic  contributions  to  the  ener,2y, 
the  total  energy  of  the  system  will  simply  be  the  sum  of  the  various 
energy  contributions.  Moreover,  a  complete  set  of  discrete  quantum 
levels  of  energy  is  aosociated  v;ith  each  conxribution.  Each  set  of 
levels  of  energy  may  be  represented  schematically  by  the  following 
diagram. 

*  The  levels  (discrete  quantum  levels)  of  energy  are  numbered 


as  indicated.  The  lo;j-est  or  ground  level  is  the  minimum 
possible,  or  allov^able,  level  of  a  molecule. 

*  The  energy  of  each  level  is  represented  by  fo /£,,--• -^fi. 

*  The  statistical  weight  (sometimes  called  desoneracy  or 
multiplicity)  is  the  number  of  individual  states  of  virtually 
the  same  energy.  It  is  given  by  the  notation  gQ,  g-j_,  , 

Si- 

Schematic  dia'^^iTa'ii  of  a  set  of  levels  of  ener,?;/  of  a  molecule 

Level  of    lie.  of  the     Energy  of  the     Statistical 
Energy      given  level    given  level       ifeight 


Si  S; 


1 


■  2  £2  S2       - 

1  £1  Si        ■       - 

0  go  Sq 

In  order  to  evaluate  the  partition  function,  2,  g.  exp(-t./KT), 
the  g.  and  f .  for  each  energy  level  have  to  be  known.  However,  it 
has  been  shown  that  only  when  the  electronic  contribution  to  the 
energy  is  significant,  will  the  discrete  values  of  6 .  and  the 
corresponding  g.  prove  their  usefulness.  The  following  are  examples 
of  values  of  the  g.  and  f .  that  are  available  in  the  literature. 


A.  The  set  of  electronic  energy  levels  of  monatomic  hydrogen 
(Eef,  2)  where  J  is  the  quantum  number  and  (cm   )  is  the 
wave  number  vjhich  is  the  enuivalent  of  an  ener:jy  unit: 


Enertjy  Level  '  g  =  (2J  +  l)       Energy  (cm"  ) 

0  2  0 

1  2  82258.907 

2  2  82258.942 

3  4  82259.272 

4  2          .   97492.198 
»  •                • 

B.  The  set  of  energy  levels  of  diatonic  hydrogen  (Ref.3) 

Energy  Level  g  =  (2J  +  l)       Energy  (cn~  ) 

0  10 

1  2  '      91698.0 

2  2  100043.0 

3  2  100062.8 


II.   Thermodynar.iic  Probability  }}  and  Maxwell-Boltzmajin 

Distribution  Law 


The  relationship  betvreen  the  thermodynanic  properties  of  a 
system  and  the  niicroscopic  structure  of  the  system,  as  proposed  by 
Planck,  is  that  the  entropy  of  a  system,  S,  may  be  expressed  by 
the  equation 

S  =  K  In  U  (3) 

in  which,  K  is  the  Boltzmann  constant  and  ii[  is  the  thermodynamic 
probability.  This  equation  is  a  postulate;  it  is,raorGover,  a 
most  important  one,  for  it  constitutes  the  bridge  that  connects  the 
maxjroscopic  thermod^/namics  vrith  the  microscopic  thermodynamics.  The 
thermodynaifiic  probability  is  defined  as  the  nuinber  of  microscopically 
different  ways  in  which  a  thermodynamic  state  can  be  realized,  or  the 


number  of  different  v/ays  in  imicli  a  statistical  state  can  be  attained, 

As  an  exa.-aple,  let  N  distinguishable  particles  be  arranged  into 
X  groups  with  H|  in  the  first  group,  il  in  the  second  group,  and  so 
on  to  N ,  in  the  last  group,  and  such  that 

.  N3_  +  Ng  + *  *^i  ^ +  U^  =  N        ■      (4) 

The  total  nvmber  of  ways  of  distributing  the  N  particles  into  x 
groups  is  Nj.  This  includes,  as  different  vrays,  the  result  of  per- 
mutations of  the  N.  within  the  i  group.  However,  permutations 
vrithin  a  group  do  not  result  in  "different"  ways  of  distributing 
the  JI  particles  into  the  x  groups.  The  number  of  such  permutations 
within  a  group  is  N.j  .  Therefore,  the  total  number  of  "different" 
wa;^'-s  of  distributing  the  N  particles  into  N  groups  is  obtained  by 
dividing  the  ITJ  ways  by  II  N.',.  This  result  is  the  thermodynamic 
probability,  If. 


II  TiA 
1* 


(5) 


After  knowing  the  thermodynamic  probability,  '/f^the  following 
statement,  shows  how  energy  level,  statistical  weight,  and  thermodyna- 
mic probability  are  related  to  the  Maxv;ill-Boltzmann  energy  distribu- 
tion. At  a  given  temperatui'o,  the  distribution  of  molecules  among 
the  possible  states  of  exiergy  is  described  by  the  Maxwell-Boltzmann 
statistics,  which  the  molecular  system  follovrs  as  proved  by  experi- 
ment. This  distribution  gives  the  average  population  of  energy  level 
£.  in  terras  of  the  pop'ilation  of  a  state  at  the  ground  level  having 
as  energy  £„. 

A  number  of  atoms  or  molecules  v?ill  be  distributed  among  energy 
levels  vrith  values  denoted  by£^,  ^  -i ,  £  p>  ^^^  ^o  °^'  Sometimes  it 
happens  that  certain  of  these  levels  are  degenerate  -  that  is,  can 
be  realized  in  more  than  one  way  or  through  more  than  one  combina^ 


tion  of  quantum  numbers.  The  nu;..ber  of  different  states  having  the 
same  energy  level  is  the  "statistical  weight"  of  that  level.  Statis- 
tical weight  is  equivalent  to  having  several  boxes  of  the  mechanical 
analogue  at  ti':e  saiae  height  in  the  gravitational  field.  The  number 
of  such  boxes  at  any  one  height,  or  the  quantum  number  corresponding 
to  any  one  level,  is  said  to  be  the  statistical  weight  of  the  energy 

levsl. 

Fnen  a  situation  of  this  kind  occurs,  it  is  necessary  that  the 
energy  be  listed  as  many  times  as  the  value  of  the  statistical  weight 
of  the  level.  To  express  the  thermodynamic  probability,  it  is  desira- 
ble to  employ  double  subscripts  for  the  iJ's.  Thus,  if  i  specifies  the 
energy  level  of  € . ,  then  N^  .  will  denote  the  number  of  particles  in 
the  jth  representation  of  that  level  if  j  talces  on  the  values  1,  2, 

,  g.,  with  g.  equal  to  the  statistical  weight.  The  thermodyna-  - 

mic  probability  of  the  system  is 

"   II  N.  I   II  II  IT.  .  I  ^  ^ 

i   1-   ;  ;   iJ' 


or 


^  =  iii?N...         ■  :.      (^^ 

5   j   10' 


and 


Si 

I   IT.  .  =  IT.  (8) 

0=1   ^'  " 

where  IT.  is  the  total  number  of  -oarticles  with  energy  €  ■  . 

With  the  foundation  shovm  above,  the  Maxwell-Boltzmann  distri- 
bution Law  can  be  derived,  vj-hich  shows  how  distinguishable  atoms  or 
molecules  are  distributed  among  their  energy  levels.  Consider  K 

particles  vrith  allovrable  energy  levels  £„,  £-,,  2  ,  etc.,  which  have 


statistical  weiirhts  g^,  g^,  s^,    etc..  The  total  energy  of  these  par- 
ticles is  a  fixed  quantity,  3.  Since  both  the  total  number  of  raole- 
culos  and  the  total  energy  are  fixed,  this  statistical  problem  in- 
volves two  conditions  of  constraint. 

By  using  the  method  of  Lagrangian  multipliers,  the  maximum 
thermodynamic  probability  can  bo  establiahod. 

Taking  the  logarithm  of  equation  (7),  there  results 


In  W  =  N  In  N  -  K  -  2  E  (l\f   In  W .  .  -  N   ) 

i  J    -'•J     -'•J    -^J 


(9) 


Moreover 


="?  ''i=?J"lO  (10) 

and  the  total  energy  constraint  is  expressed  as 

a  =  1  N.  t  .  =2.£  :■!.  .t  .  .  (11) 

Talcing  variations  of  equation  (9),    (lO),    (ll) 

i  In  ¥  =  -ZI  In  N.  .til.  .  (12) 

^B  =  l^H.    =2IXn.  .    =  0  (13) 

i  1      1      i  J     1        ij  ^    ^^ 

and  multiplying  equation  (13)  by  <x    and  (I4)byy3  ,    the  combination  of 
(12),    (13),    and  (14)   gives 

^In  vr  -o(Sl]  -^S2.  =  -It   (In  II.  .   +  •<    +iSf  .)^K.  ..        (15) 

To  render  In  ir  a  maximum,  the  coefficient  of  each  of  the  variations 
of  the  right-hand  side  equation  of  (15)  is  set  equal  to  zero.  The 


8' 


result  is  that 

In  N.  .  =  -er-fi£..  (l^) 


or 


N.  .  =  e 


-«-y8«i  (17) 


Since  N.  .  depends  only  on  i  ( throu.£;;h  €  ^)  and  not  on  j,  the 
total  number  of  particles  with  energy  t^   can  be  ifritten  as  follows 


•??. 


.-«-f«i 


"l  =  .Sj  ==,3  -  %  \,  -   %e  --^  (18) 

The  two  equations  immediately  above  are  statements  of  the  Haxwell- 
Boltzraann  energy  distribution  law. 

The  average  energy,  g  ,  of  the  particles  in  the  system  can  be 
found  as 

1  =  -^.^^" 


.«£.■ 

The  factor  e   obviously  drops  out  when  equation  (17)  is  substitued 
into  the  expression  for  the  quotient,  e/N.  It  is  evident,  therefore, 
that  A    depends  only  on  the  average  molecular  energy  and  hence  must 
be  an  intensive  veiriable. 

/3   is  identified  by  considering  its  relationship  to  the  en- 
tropy. 

3  =  K  In  If  (20) 

max 


The  entropy  for  a  system  in  equilibrium,  one  for  vrhich  the  thermo- 
dynamic  probability  is  maximum,  can  be  evaluated  by  multiplying 
equation  (9)  by  K  and  replacing  In  K.  .  by  -cr-y3£..  The  expression 


is  thus 


S  =  N  K  In  W  -  N  K  +  KiZ.N.  ,  (  «<  +j8  g.  +  l)       (2l) 


or 


S  =  K  In  N.'+  N  K«  +  K/3E  +  N  K  (22) 

Differentiation  of  (22),  while  keeping  N  and  all  of  the  . »s  cons- 
tant, gives 

■  dS  =  N  K  dc<'+  K  E  d^  +  K/3dE  (23) 


Since 


N  =  1  H  =  £  g  e  (24) 


I   -^   t 


then 

dN  =  0  =I(-do<-^.dp)  g.e"**"^^' 

=  -  £  N.  (d<x  +  £.dp  ) 

=  -II  dfi<  -  E  d/3  (25) 

Combining  equation  (23)  and  equation  (25): 

The  differentiation  of  S  is  obtained  by  holding  constant  all  in- 
dependent variables  except  3.  Prom  classical  thermodynamics, 


10 


therefore,  by  comparison  of  (26)  vfith  (2?)  "       "..  -   ; 

P^-W-  ■"        (28) 

The  Maxwell-Boltzmann  distrihution  law  equation  (l8)  can  he 
v;ritten  in  the  following  form 

_  !i  ,         ■       (29) 

N.  =  C  0  ^^ 
1 

where  C  is  the  g.e   . 


ill.  Expression  of  Thermodynamic  Properties  of  a  Gaseous 
System  in  Terms  of  Partition  Function 


1.  Internal  Energy  2 

Prom  equation  (I9)  and  (28),  the  internal  energy  S  of  a  system 
can  be  written  as 

_  ii 

Sli.   N.     S2£.p-.e  ^^ 
E  =  ■    '<   "-     "-     =         i     ^-"^^  >■ 


*  1    ti  ' 


n?€iSie 


Q 
where  Q  is  the  partition  function, 

Q  =  *  g  •  e 

Wow  taking  the  derivative  of  Q  with  respect  to  T, 


(30) 


(31) 
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dQ      1    <  «  ,.  ."  KT 


dT       2   •  i  i 

KT    ' 


^^C,o  -^  •    (32) 


Prom  oquation.T,  (30)  and  (32),  E  can  be  expressed  as 


"  "    Q     dT 


or 


MT-^^-^°^ 


2.  Entropy  S 

The  entropy  of  a  system  can  likewise  be  easily  computed.  From 

and  by  use  of  the  definition  of  Q, 

Q  =  Ne  .  (34) 

By  subsitution  of  the  internal  energy  or  given  by  equation  (33)  and 
from  equation  (34)  into  equation  (22),  and  by  use  Stirling's  approxima- 
tion, the  entropy  can  be  expressed  in  the  following  simple  closed  from: 


S  =  K  In  1^  /  +  WA<x  +   K^E  +  ilK 


=  K  (K  In  N  -  IJ)  +  M  In  Q  -  M  In  W  +  ^^(glnT  J^  ^'^ 

.  BK  (1„  Q  .  4121-  ) 

-  MC^—  (T  In  Q)    .  (35) 
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3.  Helmholtz  Free  Energy,  A  =  S  -  TS  . 


A.MT4^-BCT(in«.|^) 


=  -  MT  In  Q  (36) 

4.  Enthalpy,  H. 

=  mt(4^^  .;  (37) 

Then  the  pv  term  is  .      -■  ' 


T  (38) 

and  from  equation  (33)  aiid  (38)  the  enthalpy  H  can  be  expressed  as 

H  =  E  +  pv 

V  T** 

5.  Gihbs  Free  Energy,  G  =  H  -  T3. 

G  =  H  -  TS  •        ' 


13 


6,  Specific  Heat  at  Constant  Volume,  C  ,  and  Constant  Pressure  Cp. 

V 


Tbr^'i'^)] 


^^^^    Ll  91nT/  ^    aT  iPlnT/   J 

V  V 

^^•^UainT/         ainT  UlnT/    J 

V  V 


For  a  perfect  gas,   fi  =  pv/T  =  tlK  -j^       ,    and  Cp  =  E  +  C^  ,   hence 

V 


IV,  Calculation  of  Thermodynamic  Properties 


1.  Procedure 
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A.  Partition  t'unction  Q 

a.  Assuming  that  Q,  ,.   can  be  approximated  by  replacing  the 
summation  of  equation  (l)  by  an  integral,  upon  intergration 
the  following  equation  results,  (Hef.  4)» 


^trans  =  ^1  ^"^'  ''^ 


trans 

b.  Assuming  that  the  diatomic  molecule  is  a  rigid  rotator, 

upon  replacing  the  sum  of 

£irot 

Q  ,  =  z  g.e 

^rot   i  '^1 

the  follo.ving  equation  is  derivable.    (Ref.    4) 

^rot  =  ^2   ^'"^   ^'^^ 
^         _     8n^IKT 


where  r  is  the  interatomic  distance,  I  =yur, /w  is  the  re- 
duced mass,  m,  ♦  mp/m,  +  m_  for  diatomic  molecules,  and  o^ 
is  the  symmetry  number. 
c.  Assuming  that  the  diatomic  molecule  is  harmonic  oscillator, 

upon  neglecting  the  higher  terms  of 

ivib 


-vib   {  ^1 


the  follo-ifinj  equation  is  derivable,  (Eef.  4). 


^ib 

=  ^3  ^/'e'  -e) 

n 

-  rr          1 

^vib 

l-e-Z^ei^-e 

where /3   is  the  equilibrium  vibrational  constant,  and<o 
'  e  '      e 
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is  the  fundamental  vibrational  frequency. 

d.  The  exact  Q  -,   is  as  follows, 
elec  '  ' 

^elec  =  U   ^^i'  %^ 

_  il       ^ 

^  ^      KT      "  KT 

°-elec  =  So  +  Sie     ^  S^^  f 


©•  Q4.  4.  T  is  the  products  of  Q,    ,  Q  ^,  Q  .,  and  Q  . 
^total       -^  ^trans'  ^rot'  ^vib     ^eleo 


^total  =  ^^  ^ 


B.  Internal  Energy  S,  (Ref.  4) 


trans 

=  MT=  (41 

V 

2    3 

■(^ 

In  • 

2irrnKT 

=  -|-WI{T^  (- 

f) 

=  ^  m<:t 

=  -Y-  ^^T 

rot 

HT 

p          = 

RTx  (-—  + 

1 

-) 

^vib 

1-e 

+  In 


^) 


where  x  is  equal  to  h'^/KT,andw  is  the  fundamental 
vibration  frequency. 
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i-l 

B^K^ 

S 

e 

la 

elec 

1 

+ 

S 

tl 

KT 
e 

total 

=  IE 

s 


^0 


C.  Enthalpy  H 

total    total 

D.  Entropy  3,  (Ref.  4) 

Upon  substituting  in  all  the  numerical  universal  values  of 
constants,  such  as  h,  K,  and  Nq,  the  following  equation  is 
obtained. 

^trans  "  ^'^°^^   1°S  m  -  2.3141  +  11.4391  log  T 

39 
^rot  "  4.5757  los{^~-)       -   O.78O4  +  4.5757  log  T    >   . 

il  -fl    ■  ■• 

«     =Hln(l.ge'^^).^-^A^^/   ■  > 


1  +  g  e 


€1 
KT 


where  g  = 


^vib  =  «  ^-(^)  ^  1] 


total   ^ 


E,  Helmholtz  Free  Energy  A,  (Eef.  4) 
A  =  -  HT  In  ( -|-  +  1  j 


F.  Gibbs  Free  Snergy  G,  (Eef.  6) 


^, ''Va?;, 


G  =  -  RT  In 


N 


G.  Specific  Heat  at  Constant  Volume  C 


aE 


V  trans 


3   \3T  / 

3 


2 

3 


m 


R 


C    ,  =  R 
V  rot 


'v  vib 


R  .^  e" 


where  x  = 


KT 


^1 


'v  ele 


=  ^Ckt) 


s. 


■0 


£i 

KT 


1  + 


^0  "  KT 
e 


'v  total 


=  1  c 


H.  Specific  Heat  at  Constaiit  Pressure  C 


p  total    V  total 


2.  Exarnple  1. 

Calculation  of  the  properties  of  one  mole  of  monoatomic  hydro- 
gen H  ideal  gas  which  is  under  the  condition  of  temperature  273.16  K, 
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voliiiue  22416  c.c,  aiid  pressure  1  atm. 

For  a  monoatomic  gas,  only  trans lational  and  electronic  terms 
need  to  be  calculated.  Values  of  g^,  e^,  and  t -^  are  obtained  from 
"Atomic  Energy  Level"  N.  B.  S.  Circular  4^7  (Ref.  2). 

A.  Partition  Function  '-.'• 


•^..an.  =(-^f '  (-)^^^  ^ 


h 


.-Id 


^  3/2  3/2 


/2.x  1.38x10^  )   fi-iL^Ii^)    (22416) 


(6.26  X  10"^^)^  ^        ^  6.02  X  10' 
28 


=  1.87  X  10 


Q  +  =  1 
^rot 


^vib  =  1 


_|l      _  £2 
^elec  =  ^0  ^  Si^     ^  ^2^ 


_  £1    _  82258.9  X  1.98562  X  10  -^^ 


If  T  is  lo..,  e"  ^^  =  e~   1-38  x  10-16  ^  273.16   x  0 


Q  -,   =  2 
elec 

^total  =  11  ^i  =  (l-S'7  ^  10^^)  ^  2  =  3.74  X  10^^ 


B.  Internal  Energy  E 


^trans  =  -|-  (1-987  x  273.16)  =  8I4  cal/mole 


E   ,  =  0 
rot 


E  .,  =  0 
vib 
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elec  K       Sq 


»0 


£l        ■ 

KT 

e 

tl 

P" 

b]_ 

_ 

KT 
e 

^0 

■•  - " '.    ■■ 

E+    ,    .    =  Z  E.   i  814  cal./mole 

total         1  '  ■   •  :;   : 

C.  Enthalpy  H 

^^total  "  ^total  ■"  ^'^^^14  +  (1.987  x  273.16)  i  I356  cal./mole 

D.  Entropy  S  . -• 

^trans  "  ^'^^^^  ^^S^q^  -   2.3141  +  11.4391  Ios^qT 
=  25.7 

£         _^     KT 

3  .^^  =  H  m  (i-^-i^  e"^)  .  -I  -i^-^ ■'  ■ 

elec       ^   g-       '   KT 


1  O'  -Si- 

1  .  -^  e"  ^^^ 
Si 


a  1.38 


^Otal*  25.7  ^  1.38 

»  27.08  cal./mole  deg. 

E.  Helmohltz  Free  Energy  A 


A  =  -  ST  ln(-4-  +  1) 
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(1.987  X   273.16)  Inf  ^'^^  ""  ^^,.^     +  l]=  -  456 

^   6.02  X  10"^^     "^ 


F.  Gibbs  Free  Energy  G 
G  =  -  ET  In  -|- 

=  -  (1.987  X  273.16)  ln[-^^-24-2LJ:2^]=  _  456  cal./mole 

6.02  X  10"^^  "^ 

G.  Specific  Heat  at  Constauit  Volume  C 
C^  =  _1_H  =  -J-  (1.987) 

=  2.98  cal./mole  deg. 
H,  Speoifio  Heat  at  Constant  Pressure  C 


P 


°P=^^  =  -f^^-987) 


=  4.96  cal./mole  deg. 
Exajuple  2. 

Calculation  of  the  properties  of  one  mole  of  diatomic  hydrogen 
Hp  ideal  gas  which  is  under  the  condition  of  temperature  I5OO  K, 
volume  22416  CO.,  and  pressure  1  atrn. 

Since  Hp  is  a  diatomic  molecule,  the  trans lational,  rotational, 
vibrational,  and  electronic  terms  have  to  be  counted  in  the  calcula- 
tion. Values  of  I,  ^  ,  Wg,  E-,'  ^n'  ^1'  ^^®  obtained  from  "Spectra  of 
Diatomic  Molecule",  Herzberg  (fief.  3) 

A.  Partition  Function  Q 

2rrz  1.36x10-^^  ^^/' - -^2 

\6.62  X  10"2^)2    -'   \    5^02  X  10' 

^30 


«-3  =fyn7t1?i#-J    fi-;^3)    ^-") 


=  3.74  X  10- 
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=  2.48  X  10^^  -^ 


where  ^  is  2  and  I  is  O.46  x  10         (g^cm  ) 


«rct  ■  8-5<5 


:( ) 


^vit   °  "(    .  ,,   .'  ■         ', 


where  <•>€   is  4395.2  (era"   ) 
X  =  h*^VKT  =  4.21 


^vib  '     1  -  0.0015       '  "^'^^ 

_  €1 

*elec   "  &o  "   "  " 


KT 

^al  a.,^      ~     Sr,     *     gnS  + 


9I689   X   (1.96  X   10"-^^) 

-   1    ^     ~  K   1500 

=  1  +  e 


*»  1 

^total  =  ^I  "^  =  (3.74  X  10^°)   (8.56)   (1.02)   (1) 


= 

3.27  X  10^-^ 

B. 

Internal 

Energy  E 

trans 

—-  l^KT 

a 

-|-  2T 

=  4470 


\ot   =  i^"  =  2980 
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E,^.  =  STx  (-1 7^^ ) 

=  6270 


e  -  1 


So 

KT 
e 

Si 

^0          KT 
1  +  e 

Si 

elec      K   ~ 


E^  4.  1  =  13720  c^l./mole 
total 

C,  Enthalpy  H 

^total  "  ^total  *  ^ 

=  13720  +  (1.987  X  1500) 
=  16700  oal./mole 

D.  Entropy  S 

trans 

S  ,  =  6.53 

rot 


S  .,  =  0.2 
vib 


S  ,   =  0 
elec 

S,  .  -]_  =  42.82  cal./mole  deg. 
E.  Helmoltz  Free  Energy  A 
A  =  -  ET  In  (-4-  +  1) 

=  _  HT  m  (  ^'^T  X  10^^  ,  ,  ) 

6.03  X  lO'^-^ 


=  -  53,000  cal./mole 
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V 


F.  Gibbs  Free  Energy  G 
G  =  -  HT  In  {-^) 

=  -  53,000  cal./mole 

G,  Specifio  Heat  at  Constant  G 

C  ,     =  -4-  R  =  2.98 
V  trans    2 


C    +  =  R  =  1.98 

V  rot 

xj.  2  X 

V  vib    (e-^-  1)'^      ^  , 

il 

_f0  ~  KT 
^v  elec  ~  ^^  KT/ 


£1 

:,-:  ^0  ~  KT  :■•: 

1  +  e 

gl         •..;....■....,  ^- 

C   ,  ,  T  =  5*503  oa.l./iiiole   deg.  ',:  "  ;, 

V  total  ■'''-'  '  ^ 

H.  Specific  Heat  at  Constant  Pressure  C  '   ■,• 

C  =  C^  +  R  =  5.503  +  1.987  =  7.49  cal./mole  deg.  • 

Example  3.  Calculation  of  specific  heat  at  constant  pressure,  C  , 
and  specific  heat  at  constant  volume,  C  ,  of  Hp  as  a  function  of 
temperatures  from  6OO  K  to  6000  K. 

Calculation: 

The  total  specific  heat  at  constant  volume,  C  ,  will  be  the 


sum  of  C   ,     ,  C     ,  and  C    ., 
V  trans'   v  rot      v  vib 


C,,-,=C,     +C    ^+C 
V  total    V  trans    v  rot    v  vib 

=  ^-R  +  R  +  -^^^ 


(e--  D' 


24 


Tfhere  s  = 


cj  h 

KT 


aiid  the  total  specific  heat   at  constant  pressure,    C    ,   \iill  be 
p  total         V  total 


TABL3  I 

T  (°K-) 

Rx^o'^ 

C 

V 

C 

V 

(e--l)2 

cal/mole,des. 

cal/mole,de^. 

600 

0.0031R 

2.503IR 

4.95 

3.503IR 

6.94 

1200 

O.I489R 

2. 6489 R 

5.24 

3.6489R 

7.23   • 

1800 

0.3974R 

2.8974R 

5.72 

3.8974R 

7.70 

2400 

O.5836R 

3.O836R 

6.10 

4.O836R 

8,08 

3000 

O.7035R 

3.2035R 

6.34 

4.2035R 

8.31  ^:  :' 

3600 

O.7817R 

3.2817R 

6.50 

4.2817R 

8.48 

4200 

O.8318R 

3.33I8R 

6.59 

4.3318R 

8.57 

4800 

O.8684R 

3.3684R 

6.65 

4.3684R 

8.65    : 

5400 

O.895OR 

3.395OR 

6.70 

4. 395OR 

8.70 

6000 

0.9137H 

3.4137R 

6.75 

4.4137R 

8.80 

TABLE  II.    GAS  TABLE:    KESNAN  AM)  KAYE 


T 

^ 

C 

■D 

°K 

cal/mole,deg. 

oal/mole,deg. 

600 

5.015 

7.002 

1200 

5.415 

7.402 

1800 

6.019 

8.007 

2400 

6.472 

8.459 

3000 

6.801 

8.787 
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The  re3ulta  of  calculation  ar©  plotted  In  Fig.  1, 

and  compared  with  values  from  "Gas  Tables"  Keenan 

and  Kaye  on  the  same  figure. 
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Fig.  1  Molar  Specific  Heat  of  Hg  as  a  Function  of  Temperature, 
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3.  Discussion: 

All  the  calculations  are  based  upon  a  gaseous  molecular  model 
which  follows  exactly  the  Maxwell-Bo It zmann  distribution.  This 
gaseous  model  is  an  ideal  gas.  Hence,  the  larger  the  volume  of  the 
gas  the  smaller  the  discrepency  between  the  results  calculated  by 
this  method  and  the  propartiss  of  the  real  gas.  Moreover,  in  this 
calculation  the  following  relations  are  neglected,  1,  the  non-rigid 
rotating  property  of  a  diatomic  molecule,  2,  the  anharmonicity  of 
the  vibrating  property  of  the  molecule,  3,  the  interacted  vibrational 
rotating  property  of  the  molecule,  4,  the  electronic  vibrating 
property  of  the  molecule.  The  error  introduced  by  the  ideality  can 
be  classified  as  external  error,  while  the  errors  due  to  non-rigidity, 
the  anharmonicity,  and  the  internal  interactions  can  be  classifed 
as  internal  error,  A  more  precise  analysis  is  possible,  thus  reducing 
the  number  of  simplifications  used  but  at  the  expense  of  much  mathe- 
matical complexity. 

As  the  temperature  of  the  system  goes  higher,  dissociation  takes 
place  according  to  the  reeiction,  Hp  r   ( 1  -  y)Hp  +  2yH,  and  mon- 
atomic  hydrogen  will  be  present  in  the  system.  Moreover,  when  the 
temperature  reaches  certain  rather  higher  values,  which  are  always 
greater  than  for  dissociation,  then  ionized  hydrogen  and  free  elec- 
trons will  appear.  In  comparison  of  the  specific  heat  at  constant 
pressure,  C  ,  and  the  specific  heat  at  constcint  volume,  C  ,  calculat- 
ed by  this  statistical  method  with  the  data  tabulated  in  "Gas  Table", 
by  Keenan  and  Kaye,  the  discrepency  is  relatively  large  at  higher 
temperature.  The  reason  for  this  is  that  Hp  dissociates  appreciably 
at  about  3C00  K  and  this  effect  has  been  incorporated  into  the  value 
tabulated  in  the  "Gas  Tables", 

However,  a  more  precise  statistical  method  would  remove  this 
discrepency  and  has  the  further  advantage  of  predicting  thermo- 
dynamic properties  up  to  very  high  temperatures,  far  above  the 
experimental  ranges. 
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The  purpose  of  this  report  is  to  show  the  relation  betv^een  the 
thermodynamic  properties  of  a  gaseous  system  ajid  the  partition  func- 
tion. The  partition  function  is  defined  as  Q  =  2  g.  '^XT{-^./kt), 

j   1       1 

¥hen  the  energy  levels  i. . ,  the  statistical  weights  g. ,  and  the 
temperature  of  a  gaseous  system  are  knovm,  the  partition  function  of 
the  system  is  fixed.  In  statistical  thorraodynamios,  all  thermodynamic 
properties  of  a  gaseous  system  can  be  vrritten  in  terms  of  the  par- 
tition function,  hence,  a  thermodynamic  property  is  merely  a  function 
of  the  partition  function  from  this  point  of  view. 

In  statistical  thermodynamics,  the  concept,  S  =  K  In  Vf,  is  first 
introduced,  where  S  is  the  entropy  of  the  system,  K  is  Boltzmann 
constant,  and  \I   is  the  thermodynamic  probability.  By  maximizing  ¥, 
or  rather  In  17,  the  system  will  be  in  equilibrium.  This  is  .carried 
out  by  using  the  method  of  Lagrange  Multipliers  and  the  Stirling's 
Approximation.  Consequently,  the  entropy, 3  =  NK  ^^  (T  In  Q),  can  be 
obtained,  where  N  is  the  number  of  particles  in  the  system  and  T  is 
the  absolute  temperature  of  the  system. 

From  the  Maxwell-Boltzmann  energy  distribution  law  and  the  par- 
tition function,  the  internal  energy  can  be  expressed  as 

N£i.  N. 


t   IT. 
i   1 

and  finally,  E  =  NKT  (ilnQ/^lnT)  is  obtained.  According  to  the  defini- 
tions, all  other  thermodynamic  properties,  such  as  A,  G,  H,  C  and  C  , 
can  be  obtained  by  using  the  loiown  values  of  S  and  E. 

The  value  of  the  microscopic  statistical  method  to  evaluate 
thermodynamic  properties  of  a  gaseous  system  lies  in  two  facts;  first, 
the  thermodynamic  properties  of  a  gaseous  system  that  may  be  unobtain- 
able by  macroscopic  classical  experiment  can  be  evaluated  using  this 
method,  and  second,  even  though  the  macroscopic  experiment  can  be 
carried  out,  the  statistical  method  yields  thermodynamic  information 
with  a  higher  accuracy  and  precision  than  the  classical  method. 


